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ABSTRACT 

We construct a simple but self-consistent analytic ionization model for rapid exploration of 21cm 
power spectrum observables in redshift space. It is fully described by the average ionization fraction 
Xe{z) and HII patch size R{z) and has the flexibility to accommodate various reionization scenarios. 
The model associates ionization regions with dark matter halos of the number density required to 
recover Xe and treats redshift space distortions self-consistently with the virial velocity of such halos. 
Based on this model, we study the line-of-sight structures in the brightness fluctuations since they are 
the most immune to foreground contamination. We explore the degeneracy between the HII patch size 
and nonlinear redshift space distortion in the one dimensional power spectrum. We also discuss the 
limitations experimental frequency and angular resolutions place on their distinguishability. Angular 
resolution dilutes even the radial signal and will be a serious limitation for resolving small bubbles 
before the end of reionization. Nonlinear redshift space distortions suggest that a resolution of order 1 
- 10" and a frequency resolution of lOkHz will ultimately be desirable to extract the full information 
in the radial field aX z ^ 10. First generation instruments such as LOFAR and MWA can potentially 
measure radial HII patches of a few comoving Mpc and larger at the end of reionization and are 
unlikely to be affected by nonlinear redshift space distortions. 
Subject headings: cosmology: theory — diffuse radiation — methods: analytical 
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1. INTRODUCTION 

The epoch of reionization is perhaps one of the 
"darkest" place in modern cosmology now adays. On 
one hand, Gunn-Peterson constraints (Gtmn & Peterson 



of lin ear velocity flows o n 21cm power spectra fBarkanaJ 
|2005t IBarkana fc Loebl l2005ald: iBharadwai fc All .20051 



196S) from quasars (e.g.. iBecker et a/.. ,2001; .Fan et 



2002(1 tells us that the Universe was mostly ionized by 
z ~ 6. On the other hand, the optical dep th to Thom- 
son scattering r "^ 0-17 f rom WMAP data l|Kogut et al\ 
1200."^ ISoergel et oiJl200.1 1 implies that reionization may 
have begun as early as redshift z ~ 17. While these data 
suggest reionization was probably an extended and com- 
plicated process, we still have little information about 
what exactly happened during this period. 

Among various approaches to explore the reionization 
epoch, 21cm radiation from neutral hydr ogen is in prin- 
ciple o n of the most p o werfu l t ools (e.g. iHogan fc Reef 
12Z| IScott fc Ree^ IMadau. Meiksen fcR.ees 

19971) and has recently received much attention 
in preparatio n for a new generat i on o f experi- 
ments (e.ff., IBarkana fc Loebl l2005a| ICalim 2005; 



Chen fc Miralda-Escudel 12004 Iciardi fc M^ 



20031: iFurlanetto. So kasian. fc Hernauis 
Furiane^g^^Zaldamaga, fc Hernquisj ^ 
Furlanetto. McOuinn, fc Hcrna uisTl 120051: IMorales 
2005; Morales fc Hewitt 20031: iPen. Wu. fc Pete~ 



20041 IS^s. Coorav. fc Knoxl l2005t IWvithe fc Loeb 
2004bl: IZaldarriaga. Furlanetto. fc Hernauis^ I2004j) . 
Although contamination from foregrounds will likely 
dominate the total signal, their smoothness in frequency 
should allow the measurement of radial structures in the 
brightness fluctuations. 

Radial structures however are measured in redshift 
space and are distorted by the peculiar velocity of the 
neutral gas. Recent studies have addressed the impact 



iDesiacaues fc Nusseill2004|l . 

In this paper, we provide a complete treatment includ- 
ing the non-linear redshift space distortions or "Fingers 
of God" that blur out radial structures along the line 
of sight. We show how these effects impact 21cm power 
spectra, how they manifest themselves with different ex- 
perimental setups, and how they interact with intrinsic 
reionization properties such as the size of the HII regions. 
We show that nonlinear redshift space distortions will ul- 
timately limit our ability to study small HII structures 
during reionization. 

The outline of the paper is as follows. In Section |21 
we construct a simple but self-consistent analytic ion- 
ization model based on an association with dark mat- 
te r halos. Our model is similar in cons truction to that 
of IFurlanetto. Zaldarriaga. fc HernauistI ll2f 



Mai) but al- 
lows for a rapid exploration of the possible parameter 
space of alternate models. We discuss a generalization of 
our model that includes a HII region size distribution in 
the Appendix. In Section 13 we employ the association 
with dark matter halos to study redshift space distortions 
based on their bias and virial velocity. In Sectional we 
study the impact of experimental angular and frequency 
resolution on resolving the features induced by the HII 
regions and redshift space distortions. In Section [SJ we 
discuss the impact of these considerations on first gener- 
ation experiments such as MWA^ and LOFAR^ as well as 
the design of an experiment that can recover essentially 
all of the information in the radial temperature field. We 
summarize our results in Sectional 

1 http://web.haystack.mit.edu/MWA/MWA.html 
^ http://www.lofar.org 
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2. REAL SPACE 21CM POWER SPECTRUM 

2.1. Brightness temperature 

The redshifted brightness temperature fluctuation in 
the 21cni hne can be viewed as a three dimensional field 
where the position is specified by the angular coordi- 
nates on a small patch of the sky and the observation 
frequency v. In the absence of redshift space distortions, 
the mapping to comoving coordinates 

r::^ Dj_{z)e + D{z)er (1) 

involves the redshift z = 1^21/1^ ^1, where 1^21 = 
1420.4MHz, the comoving transverse (or angular diam- 
eter) distance D±{z), and the comoving distance along 
the radial direction eV, D{z). In a flat universe, which 
we assume throughout, D± = D. 

The amplitude of the brightness temperature fluctua- 
tion is proportional to the neutral hydrogen density fluc- 
tuation dm = Snm/nii 



5T{f)^F^/-\z)5m{i^. 



(2) 



where is the mean hydrogen density. Under the as- 
sumption that the spin temperature is much larger than 
the CMB temp erature, the proportionality coefficient 
l|HoEfan fc R,eesllT979j) 



F^l\z)- 



3 c^hA2in]i 1 
16 kvl-^H{z) 1 -I- z 



i0.023K 



i-y„fib^^ fi + z 0.15 



0.75 0.02 V 10 ^mh"^ 



1/2 



(3) 



Here A21 is the Einstein coefficient for spontaneous emis- 
sion and Yp is the primordial helium mass fraction. 

Although the mapping involves integrals over the ex- 
pansion rate H{z), which is currently uncertain at low 
redshifts due to the dark energy, at high redshifts it is 
already well determined by CMB observations. Wc will 
consider a small range of observing frequencies around a 
central value z ^ 1. Differences in observing frequency 
are then simply proportional to differences in radial dis- 
tance 



5D 



5z 
H{z) 



17Mpc 



1 + z 0.15 
10 VL^K^ 



1/2 



5v 



IMHz 



(4) 



We set c = 1 where no confusion will arise. The trans- 
verse distance in a flat universe can be scaled forward 
from the CMB determined distance to recombination 
-D* = D±{z^,) to avoid ambiguities due to the dark en- 
ergy 



dz 
H{z) 

4.90Gpc 



1 - 



0.15 10 

rim/l2 1 + ^ 

1 + 



1/2 



1 + 2 



1 + 1 + 



1/2' 



(5) 



where the redshift of matter radiation equality is given 

by 



l + z, 



cq 



3600 



0.15 



(6) 



Notice there is no dependence on and no sensitivity to 
the dark energy. 

Likewise the conversion factor F from neutral hy- 
drogen density to brightness fluctuation only adds a 
sensitivity to the baryon density Slb/i^. Constraints 
on these quantities from the flrst year WMAP data 
ijSpergel et a/.ll200^ . require n^K^ = 0.14 ± 0.2 and 
rJb/i^ = 0.024 ±0.001 which corresponds to = 13.7 ± 
0.5 Gpc, z* = 1088^2- For definiteness and in accor- 
dance with large scale structure data sets l|Seliak et a,l\ 
|2005 iTtegiiiark aL■2004^1 ■ we will illustrate our results 
in a fiducial ACDM cosmology with fim = 1 — J^a — 0.29, 
r2^/i2 = 0.0244, rii„/i^ = 0.15. We assume that the initial 
spectrum of curvature perturbations C takes the form 



27r2 



P, 



k 



0.05Mpc" 



(7) 



with 5^ — 4.8 X 10 ^ and n = 0.99, corresponding to cts = 
0.9. This requires £>* = 13.4Gpc, = 1088 and the 
total Thomson optical depth during reionization to be 
T — 0.103 or a lower limit on the redshift of the beginning 
of reionization of z = 12. In the rest of this section, we 
shall treat the observations as if they were made directly 
in 3D physical space on the neutral hydrogen density 
field. 

2.2. Neutral hydrogen halo model 

The neutral hydrogen fluctuations can be separated 
into contributions from the ionization fraction x{r) — 
Xe + 5x(f) and gas density fluctuations 5{r) 



5m{r) = [1 - - 5x{r)\5{r) - 5x{r) 



(8) 



The two point correlation function of 
the neutral hydrogen field then becomes 
ijFurlanetto. Zaldarriaea. fc Hernauistll2004j) 

^.SmiSmA'r) = S.SxSx{r) - 2(1 - Xe)£,sxs{r) 

+ {^-Xef£,ss{r)+S.SxSSxs{r). (9) 

Here we have employed the notation 

Cab(|ri -f2|) = (a(fi)6(f2)) (10) 

for the two point correlation function between field a and 
b. Analogously we define the power spectra as the Fourier 
transform of the correlation functions 



Pab{k) 



d^k 

(2^ 



(11) 



It is also useful to define the ID line of sight projected 
power spectrum since the radial structure of the field is 
critical for removing foreground contaminants. In gen- 
eral it is defined as 



d^k 



^Pab{k) 



(12) 



where fc|| and arc the components of k parallel and 
perpendicular to the line of sight. Finally we will occa- 
sionally use the shorthand notation Pa = Paa where no 
confusion will arise. 
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The final term in equation (|5J involves a product 
of density and ionization fluctuations. Though for- 
mally second order for linear fluctuations, we shall 
see that for density perturbation on scales smaller 
than a typical ionized region, this term must be in- 
cluded to construct a physical model of the neutral 
hydrogen field because of contributions of the form 
(jFurlanetto. Zaldarriaea. fc Hernauist)l2004|) 



^SxSSxs{r) ~ ^SxSx{r)^s5(,r) . 



(13) 



We model the underlying density and ionization cor- 
relation functions with a halo model of ionization bub- 
bles . Our implementation repres ents a si mplific ation of 
the iFurlanetto. Zaldarriaga. fc H ernauistI l)2004|) model 
in that we take a single characteristic bubble size at each 
redshift. On the other hand, we leave the evolution of the 
bubble size and mean ionization arbitrary to ensure flexi- 
bility. For example, in our model reionization can end by 
the percolation of many small bubbles or the growth of 
rare bubbles until overlap. More importantly our model 
captures, in a simple and fully analytic form, the required 
physical scaling of the correlati ons as a function of mean 
ionization and bubble size (cf. ISantos et aiMOO^ . 




Fig. 1. — Mass threshold and virial temperature (see equa- 
tion lijYI ) as a function of the mean ionization Xe for several choices 
of the bubble radius ij at 2 = 12 and 2 = 8. 



By associating the seeds of the ionization bubbles with 
massive dark matter halos of that number density, we 
can model the statistics of the bubbles. Specifically we 
set a mass threshold such that 



rib = 

where dn / d In M 
(^Sheth fc To rmcn 

duh 



dM drih 
M d\nM 



(16) 



Mth 

is the halo mass function 



d\nM 



PmO r, ^ 

M ■'^'^'dlnM 



Here is the matter density today, i' 
and 



(17) 

Sc/cTliniM;z) 



A\l-aiy^[l + [aiy'yP] exp[-aj.72] , (18) 

and (Tiin is the rms of the linear density field smoothed 
by a tophat of a radius that encloses a mass M. A 
choice of 5c = 1.69, a ~ 0.75, p = 0.3, and A such that 
J di'f{i') = 1 fits the results of simulations well. Note 
that our model can accommodate collective effects where 
smaller halos associated with the seed halo can make the 
bubble ra dius much larger than the virial radius of the 
seed halo (! Barkana fc Loebi l2004'). In Figure ^ we show 
the mass threshold Mth as a function of Xe for several 
choices of the bubble size 7? and redshift. 




- (b) z = 8 



0.01 



In our model, the correlation functions are described by 
two quantities: the typical radius of ionization bubbles 
R{z) and the mean ionization fraction Xe{z). We assume 
that the probability that a given point in space is ionized 
is determined by a Poisson process such that 



{x(f)) = 1 — e 



-n(,(f)Vi, 



(14) 



where rif, is the number density of bubbles and Vj, = 
47ri?^/3 is the volume of the bubbles. Brackets here de- 
note averaging over the Poisson process. The mean num- 
ber density is related to the mean ionization by 



rib 



(15) 



Fig. 2. — Bubble bias as a function of the mean ionization Xe for 
several choices of the bubble radius ij at 2 = 12 and 2 = 8. The 
bubble bias increases with R at fixed Xe since the bubbles become 
rare and highly correlated. 



2.3. Two bubble correlations 

For scales where r ^ R, the two point functions are 
dominated by the correlations between two separate bub- 
bles as in the two halo regime of the halo model. These 
correlations are induced by the enhanced probability of 
bubble formation in over dense regions through equa- 
tion ifT^ . 
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Fluctuations in the bubble number density will then 
follow the density fluctuations of the dark matter through 
the bubble bias 



rib JMth M ' d\nM 



(19) 



where the halo bias is given by llSheth fc TormenllT999l) 

Specifically, we model 

nb{r) ^ fibil + bSw) , (21) 

where dw is the density fluctuation field smoothed by a 
top hat window of radius R, 



(22) 



where WR{r) = VJ,~^ for r < i? and otherwise. In 
Figure 13 we show the bubble bias as a function of Xe for 
several values of R and z. 




k (Mpc-i) 

Fig. 3. — Neutral hydrogen power spectrum at 2 = 8 with 
Xe = 0.8 and R = 0.3Mpc//i broken up into components from the 
two bubble and one bubble contributions. The one bubble contri- 
butions are further shown with th e fu ll convolution factor (solid) 
and the approximation of equation 1351 (long dashed). Oscillations 
at high k are an artifact of taking a single bubble scale with a top 
hat profile and should be removed by averaging (see Appendix) or 
projection. For reference the density power spectrum is also shown 
(short dashed). 

Expanding equation H14|l in the density fluctuations 
and taking the Fourier transforms yields 

P!^sAk) = [{l - xe)ln{l - x,)bWHik)fPsik) , 
5.5ik) = -{l - Xe) ln(l - Xe)bWR{k)Ps{k) , (23) 



p2b 



where Wii{k) is the Fourier transform of the top hat win- 
dow and the superscript 2b denotes two bubble contribu- 
tions. Here we have assumed that the gas density fluc- 
tuations trace the dark matter fluctuations in the two 
bubble regime. Furthermore, in this regime Pg^g g^g is 
second order in the density fluctuation and we therefore 
set it to zero. The power spectrum of the two bubble 
contributions then becomes 



where the total effective bias includes both density and 
ionization fluctuations 



6off = ln(l - Xe)bWR + 1. 



(25) 



Note that this expression has the proper limiting behav- 
ior. As the mean ionization Xe 0, the ionization bub- 
bles disappear, bes 1 and P|^j ^ Ps- As — > 1 all 
of the neutral hydrogen disappears and P^^^ 0. Note 
that this behavior is independent of R. Full reioniza- 
tion can occur either by the growth of small bubbles into 
large ones or the creation of a high number density of 
overlapping small bubbles. Finally the effective bias bcB 
of the neutral hydrogen fleld is typically negative in the 
two bubble regime given a model such as ours where ion- 
ization occurs in over dense regions. An example of the 
two bubble power spectrum is shown in Figure 01 

2.4. One bubble correlations 

For scales where r <C i?, the correlation functions are 
dominated by the presence or absence of a single bub- 
ble. Two points separated by r <C i? are either a part 
of an ionization bubble or not and hence the ionization 
probability must converge to Xe- Two points separated 
hy r ^ R have their ionization probability converge to 
x'^. Hence the ful l correlation function must interpolate 
between the two (jGruzinov 



(x(fi)x(f2)) =xl + [xe - xl)f{r/R) . 
The interpolation function f(r/R) must satisfy 

1, 



lim firlR)- 
lim fir/R)- 



■0. 



Therefore the ionization fluctuations become 
^lls. = i^e - xl)f{r/R) . 



(26) 



(27) 



(28) 



Since inside a bubble the medium is fully ionized, the 
density ionization cross correlation is negligible in this 
regime £,g^g ~ 0. Note however that ^.g^gg^g is of the 
same order as ^gg and can not be neglected. 

To determine f{r/R) for our assumed top- hat ioniza- 
tion bubbles, we consider the case where bubble overlap 
is negligible. By analogy t o the one ha lo regime of the fa- 
miliar halo model (see e.g. lCoorav fc Sheth,..2001f) where 
the halo density proflle is replaced with the ionization 
profile, we can immediately obtain f(r/R) as the convo- 
lution of two top-hat windows Vf,Wfl(r) 



fir/R) =l-|(i^) + ^(r)^ r<2R, 
= 0, r>2R. 



(29) 



We have verified through the procedure outlined in 
iFurlanetto. Zaldarriaea. fc Hernaui'g3 1)2004|) that this re- 
mains a reasonable approximation for the overlap regime 
as it simply represents an interpolation between the two 
physical limits imposed by equation l|2()|l through the 
scale R. 

Summing the contributions in equation ||5J) and tak- 
ing the the Fourier transform, we obtain the one bubble 
contributions to the power spectrum 



Pg^M^il-XeYbiMk) 



(24) 



Ps'ik) = ixe~xi) 



VbWUk) 



Psik) 



(30) 
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where 



(2^ 



(31) 



This construction guarantees the proper hmiting behav- 
ior of P}'^ as Xe ^ and ^ 1. Note that we 

OHI '^^ _ 

have arbitrarily associated the density fluctuation term 
with the two bubble contributions so that Pj^** = 0. 

The first term in equation Ij^UI) represents the shot 
noise contributions of the bubbles. In the limit where 
bubbles overlap negligibly Xe ~ nbVb 3> Xg, the first term 
becomes x\jnb as expected. The second term comes from 
Ci5s(52;('")^(5(5(^) and involves a convolution of power spec- 
tra. Note that 

^Wlik^) ^ P,(k) . (32) 

On small scales the total contribution involving the den- 
sity fluctuations directly become 

(1 - a;,)2p5(fc) + (xe - xl)P,{k) = (1 - x,)P,{k) . (33) 

Again, the reason is that for r <C -R the region is either 
ionized or neutral so that density fluctuations contribute 
to neutral hydrogen fluctuations with only one factor of 
(1 — a;^) representing the volume fraction that is neutral. 
Without this term, the neutral hydrogen power spectrum 
would be unphysical in the one bubble regime. In the 
opposite limit, 

Xm, Ps(k)^Vb I j7^Wl,{k,)Ps{ki) 



(2^)3 



= Vbal, (34) 

where tr^ is the variance of the density field smoothed by 
a tophat on a scale R. We find that a simple interpolation 
between the two regimes 



Psik) 



Psik)Vbaj 



(35) 



[(P,(fc))2 + (Ha|)2]i/2 

suffices for obtaining power spectra to 10-20% accuracy 
(see Figure El. 

2.5. Fiducial model 

The total neutral hydrogen power spectrum is simply 
the sum of the one bubble and two bubble terms in equa- 
tion (|24|l and equation l|3()|l . The 21 cm brightness tem- 
perature power spectrum in real space is then 

PsT{k)^F[P^^^{k)+P^^^{k)] , (36) 

and is parameterized by the bubble radius R and mean 
ionization Xe at a give n redshift (see Figure P| . 

Like the model of ISantos. Coorav. fc Knoxl l|2005j) . 
our model is analytic and allows a rapid exploration 
of parameter space. Furthermore l ike the model of 
iFurlanetto. Zaldarriaga. fc Hernaui'sil l)2004aD . it has the 
right qualitative behavior as a function of Xe and on 
scales smaller than the bubble radius. 

Despite the flexibility of our model, not all choices of 
Xe{z) and R{z) lead to physically plausible reionization 
scenarios. Reionization can only proceed if the seed halo 
is sufficiently large for cooling and fragmentation to oc- 
cur. For atomic line cooling, the virial temperature of 
the halo 



lO^K 



= 1.1 



0.15 



M 



10 / VlO^Mp 



(37) 




Fig. 4. — The condition that bubbles are seeded in halos of 
Tvir > 10* K sets a bubble number density and hence a bubble 
radius R given the mean ionization Xe. 



must be greater than ~ 10 K; for molecular hydro- 
gen cooling, M > 10^h~^Ma^ but the reionization 



efficiency is lower llBarkana. Haiman. fc Ostrikei l2001i 



lYoshida et all l25MISomerville. Bullock fc Livid l2003(K 

Since we are mainly interested in redshifts z 10 it is 
reasonable to choose the typical bubble radius R to corre- 
spond that required by Tvir(Mth) — 10*K and Xg. This is 
displayed in Fig.01 The fiducial model is then described 
by a single function Xe{z). 



3. REDSHIFT SPACE POWER SPECTRUM 

3.1. Redshift Space Distortions 

Redshift space distortions due to the peculiar velocity 
of the neutral hydrogen change the mapping between ob- 
servation frequency and radial distance. For the linear 
velocity field, convergences in the velocity field are as- 
sociated with over dense regions and hen ce enhan ce the 
apparent clustering of the density field l)Kaiseii ri987'l. 
The only difference between the neutral hydrogen field 
and the familiar density field example is that the neutral 
hydrogen can be an anti-biased tracer of the density field 
in the two bubble regime (see Section [^J. In this case, 
redshift space distortions suppress rather than enhance 
the power. In general the redshift space power spectrum 
P'^ of a biased tracer X of the linear power spectrum 
becomes 



P^{k) = L\hJk. 



bx)b\Ps{k) 



where fcii is the line of sight component of k and 



L{p,b) = l 



dluDc 
din a b 



(38) 



(39) 



Here Dg is the linear growth rate of the density field. For 
the redshifts we consider, Dq ex. a and dlnZ^c/lna ~ 1. 

On small scales, peculiar velocities are associated with 
the virialized motion of the gas leading to a suppression 
of radial power known as the Fingers of God (FoG ). We 
adopt a halo model for the FoG suppression (^Shethill996t 
lWhitell200i|) . The ID velocity dispersion of a halo of 
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0.1 1 

k (Mpc-i) 

Fig. 5. — Comparison of the Fingers of God fitting functions Gp 
(pairwise; exponential) and G(, (bubble; Gaussian) and the halo 
model prediction. The functions are normalized at G = 0.1 such 
that at 2: = 8, CTp = 29km/s and for bubbles with R = 0.3 Mpc/h 
and Xe = 0.8, o-;, = 22km/s. 



mass M 



1/3 



(40) 



where Avir — IStt^ is the virial over-density at high red- 
shift. This velocity dispersion translates into smoothing 
across a radial spatial scale of 

Sz 1 



H{z) aH{z) 

'flmh'^ 1 + Z 



:0.245Mpc 



0.15 10 



-1/2 



30km/s 



.(41) 



In Fourier space, the convolution becomes a multiplica- 
tion by a suppression factor of 



Wa{k\\) ~ exp 



m/aH)^al{M) 



(42) 



Halos of a range of masses contribute to the density and 
ionization power spectra and so the total suppression de- 
pends on a weighted average of these factors. 

As indicated above, we use a constant velocity disper- 
sion for each halo because the velocity dispersion profile 
has little impact on redshift space distortions. The differ- 
ence between modeling velocities through an isothermal 
assumption or through the Jeans equation in the true 
dark matter profile are small. Outside the halo scale, all 
that matters is the overall mean velocity dispersion of the 
halo. Inside the halo scale, there is only a minor differ- 
ence, the effect of taking a re alistic ah (r) is a rescaling of 
a con stant an within ±10% ({Tinker. Weinberg fc Zhenel 
l2()05l) . 

For the density fiuctuations, we separate the one and 
two halo contributions as usual 



p2hs 



plh s 



(43) 



where 



p,"-(fc>Fr(fc) 



(44) 



p2h; 



(fc): 



dM 



M 

PmQ 



dUh 



■y^k;M)WUkii).. 



dlnAf 

where y is the Fourier transform of the gas density profile 
normalized such that y{0) — 1. We take the gas density 
to trace the dark matter density in the translinear regime 
considered he re. The dark matter density pro file can be 
described by l|Navarro. Frenk. fc Wh'it^l2004D 

Pir) oc — j-^ r2 , (45) 

rc/rvir(l + rc/rvh-r 

where the virial radius encloses the halo mass M 
at the virial over-d ensity Avir and the concentration 
l|Bullock et aZ.ll200H) 

c(M) = V" " , (46) 



1 



where is defined by (Tiin(M*; z = 0) = 1.698. 
Let us define the suppression factor as 



Gp(fc) 



Psik) 



(47) 



Psik) ' 

where the real space power spectrum follows equa- 
tion H44() with Wr, 0. Th e suppression can be well 
approximated by l)Shethlll996fl 

""^^'^ " i + (fcii/Uv^ ^''^ 

or the functional form taken by an exponential distribu- 
tion of pairwise velocities with dispersion ap. To extract 
CTp, we match the numerical calculation to the fitting form 
at Gp = 0.1. An example of the accuracy of the approx- 
imation is shown in Figure El 

We follow the same procedure with FoG redshift space 
distortions of the one and two bubble contributions to 
the ionization fluctuations. The only difference is that in 
the integrals over masses in equation H44|l the lower limit 
is set to Afth and the density profile y(fc; M) is replaced 
by the bubble profile Wii{k). Again we can define the 
suppression factor as the ratio of redshift to real space 
power spectra. Since most of the contributions come 
from bubbles near the mass threshold, the suppression 
factor is better fit to a Gaussian 



Gb{k) — exp 



l{k\\/aHfai 



(49) 



We find ail by matching this to the numerical results at 
Gb — 0.1 and show an example in Figure El 

Combining the linear redshift space distortions with 
the ionization and density fluctuation FoG distortions 
gives the final form for the redshift space power spectrum 

PLik) = (1 - x,fhlsL\^i, ks)Gp{k)P{k) (50) 



+ {xe - xl)Gp{k) 



L'{ii,l)P{k) + VbWl 



where the effective bias is 

ks = \n{l-Xe)bWR+l (51) 

and we have rescaled the window function to absorb the 
differences in the FoG suppression for density and ion- 
ization fluctuations 

Gfc(fc), 



Wl{k). 



Gp{k) 



(52) 
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Note that the full redshift space power spectrum of the 
temperature fluctuations 

PiAk) = FPi^^ik) (53) 

is still defined by the cosmology and two parameters 
R, Xe- With a restriction to halos that can cool by atomic 
lines in the fiducial model (see Subsection 12. 5|l . this re- 
duces to a single parameter Xe per redshift. Conversely, 
our functional form can be used to model an even wider 
range of physical conditions if {xe, R, cTp, crt,} are all taken 
as free parameters of the model and the bubble profile 
Wn generahzed. Note that if W^j is adjusted, VJ, should 
be consistently modified. Finally, as discussed in the Ap- 
pendix, to account for a distribution in bubble sizes one 
can replace Wr with (Wb) in two bubble terms and 
with (W^) in one bubble terms which would also elimi- 
nate oscillatory artifacts in the power spectrum. 

3.2. ID power spectrum 

Foreground contamination is perhaps the largest 
problem facing all 21cm experiments. Since fore- 
grounds tend to be smooth in frequency, it has long 
been realized that one of the most robust observ- 
ables will be the radial structure i n the temper- 
ature field (Hogan&Rces 1979; Scott fc Ree'3 119901 
Gnedin & Shaver 2004; Morales. Bowman. & Hewitt J 
2005.: Oh fc Mack .2003; .Wang et al. 2005). We will 
therefore focus on the ID power spectrum for the rest of 
the paper. This will also enable us to understand intu- 
itively the impact of redshift space distortions on future 
experiments. The modeling of the 3D power spectrum 
in the previous sections, however, can also be applied to 
interpret foreground-cleaned data where the transverse 
structures can be extracted. 

Employing the general relationship between 3D and 
ID power spectra from equation (|12|l . we show several 
examples of the real and redshift space ID power spec- 
trum in the 21cm temperature fluctuations in Figure El 
We plot the variance of the temperature fleld contributed 
per logarithmic interval in k or frequency 

which then can be simply interpreted as the square of 
the amplitude of the signal. All curves have the same 
parameters except the HII bubble size R. Note that by 
flxing Xe and varying R, we are relaxing the constraint 
that bubbles are seeded by atomic cooling halos in all 
but the R = 0.3 h~^Mpc curves (see Figure0J. 

In real space, the bubble size R places a feature in 
the ID power spectrum associated with the scale of ion- 
ization fluctuations. For a large bubble radius i? > 1 
h^^Mpc and high mean ionization, this feature should 
be distinguishable from the smoother density fluctua- 
tions even in projection. In practice, this feature may 
appear as a smooth bend in the power spectrum if the 
distribution of bubble sizes is wide (see Appendix). Note 
that the signal is expected to continue to rise to smaller 
scales due to density fluctuations outside of the bubbles. 

In redshift space, peculiar velocities impose a second 
scale that cuts off the ID power spectrum regardless of 
the bubble size. We define the cut off scale corresponding 
to fjp and tJb as AD^r according to equation 141|l . For the 



density fluctuations ap controls the FoG effect, ap > <Th. 
All radial power on scales smaller than this cutoff scale 
will be erased. 




0.01 0.1 1 10 100 

k [1/Mpc] 



Fig. 6. — Real (dashed) and redshift (soUd) space ID power 
spectra for different fiducial Hlf patch sizes with perfect angular 
and frequency resolutions. All parameters are the same except Hit 
patch size R. The curves for R = lO/i'^Mpc and R = 0.03ft,"^Mpc 
have been multiplied and divided by 10 respectively for clarity. 
Note the curves represent the variance of the temperature field in 
mK^ contributed by a logarithmic interval in k. 

The existence of two scales in the power spectrum 
means that the detection of a single feature should not 
automatically be associated with the bubble scale. When 
R > ADfT, HII structures always suppress the power 
spectrum at lower k than the nonlinear redshift space 
distortion does, so that we can distinguish the HII bub- 
ble size signature from the power spectrum. With suffi- 
cient radial or frequency resolution we will also be able 
to observe that FoG distortions lead to a reduction in 
power that continues to increases with fc. 

The observation of a feature followed by a cutoff in 
the ID power spectrum will ensure that the feature is 
from HII patches. By measuring its evolution in redshift 
we obtain a statistical measure of the evolution of the 
typical size of the ionized region. However in the fiducial 
atomic cooling model with R = 0.3 h~^Mpc only the 
FoG feature is prominent. In scenarios where Xe is also 
adjusted, the prominence of the bubble feature can be 
adjusted but a confusion between the two scales can still 
occur. 

If only a single feature is measured and barely resolved 
then the ambiguity between FoG distortions and bubble 
size must be broken by measuring the transverse power 
spectrum. Fortunately, the cosmological redshift space 
distortion or Alcock-Pacynski effect can be co nsidered 
known and does not introduce an ambiguity fc.f. lNusser I 
120051) . Foreground contamination in the transverse di- 
mensions however will have to be controlled. 

4. OBSERVED POWER SPECTRUM 

4.1. Modeling experimental resolutions 

The observed 21cm power spectrum will also be dis- 
torted by instrumental effects from finite angular resolu- 
tion A6 and frequency resolution Av. These modify the 
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3D power spectrum as 

P^^^k, M, A^.) - PlT{k)Wlo{k^)Wl,{^) , (55) 

and consequently the projected ID power spectrum 
through equation (|12|l . We will take the smoothing due 
to frequency resolution to be a Gaussian of FWHM lS.v 
so that 

W-i,(fc||)=e-'^l*^ , (56) 

where 5D is given by equation with 5v ~ Ai'/\/81n2 
to convert the FWHM to a Gaussian width. 

For the angular resolution, a single dish experiment 
with a Gaussian beam of FWHM resolution A6' would 
produce a window function 



-fci_D5_(Ae/\/8ln2)^ 



(57) 



where the angular diameter distance is given by equa- 
tion ISJ. In our fiducial cosmology, this factor is approx- 
imately 



D±Ae 



V8hr2 



1.65Mpc 



1 - 0.33 



10 

l + z 



1/2 



A9 



(58) 

where we have also converted the angular units from ra- 
dians to arcminutes on the rhs for convenience. 

All planned 21cm experiments are however interfer- 
ometers. In this case the window function is given by a 
sharp cutoff in kj_ space defined by the longest baseline. 
An interferometer array with a maximum baseline of L 
measures the transverse power spectrum out to 



^cut — 



27rL 

Xd7 



(59) 



where A is the observation wavelength. With uniform 
coverage of the baselines, the window function can be 
approximated as 



wUk^) 



r 1 

I fc^ 



!i ^CUt J 

> ^cut ■ 



(60) 



Since the net effect of this sharp cut on the ID power 
spectrum is similar to a Gaussian beam of FWHM 

A 



A0 = \/8 In 2 
= V81n2 



2nL 
1 



kr^t.D I 



(61) 



we will use the Gaussian window for illustrative purposes 
in the following sections. 
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Fig. 7. — Ratio of redshift space power spectra with finite fre- 
quency resolutions Pg^°^^{k; A9, Au) to the one with perfect fre- 
quency resolution Pg^°^''{k; A8,0), for a model with parameters 
z = S, Xe = 0.8 and R = O.Sft-lMpc. 
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Fig. 8. — Ratio of redshift space power spectra with finite angu- 
lar resolutions P^^°^^ {k; AO , Au) to the one with perfect angular 
resolution P^J^°^^{k\0, Av), for a model with parameters 2 = 8, 
Xe = 0.8 and R = O.S/i-^Mpc. 



4.2. Observed ID power spectra 

The impact of frequency resolution on the ID power 
spectrum is simple in that equations H12I) and H55|l imply 

P^^°''%k;A9,A,.)^WLik) J ^wUk^) 

xP!A\/k^ + kl), (62) 

so that the relative effect of frequency resolution is inde- 
pendent of the neutral hydrogen model. In Figure [7| 
we plot W^^ or the ratio of 21cm redshift space ID 
power spectra for different channel widths Av relative 
to Ai> — 0. Comparing this figure to the features from 



the ionization bubbles and the FoG distortions one might 
naively suppose that the criteria for resolving these de- 
pend on Ai/ alone. 

However for projected ID power spectra, the angular 
resolution plays an important role. Only modes whose 
three dimensional wavevector k have transverse compo- 
nents that can be resolved contribute to the radial fluctu- 
ations in equation H62|l . Thus a finite angular resolution 
degrades the ID signal on all scales and additionally can 
lead to a new feature in the ID power spectrum. Fig- 
ure |S1 plots the ratio of redshift space power spectra with 
different finite effective beam sizes Ad to the one with 
A^ — 0. Note that the suppression effect is independent 
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of IS.V but does depend on the model 3D power spectrum. 
Here we have assumed the fiducial R = 0.3 /i~^Mpc, 
= 0.8 model at z = 8. The feature in Figure |H1 near 
1 Mpc"^ is due to the assumed bubble scale. Even with 
perfect frequency resolution, the angular resolution can 
mask features from the bubble scale and FoG effects. In 
the context of foreground removal, Figure |H1 implies that 
for > 30" the experimental angular resolution will 
place a more important limitation on cleaning than FoG 
redshift space distortions. 

Note that the total suppression is simply the product 
of the two curves in Figure and |H1 The optimal fre- 
quency resolution to extract most of the power in the ra- 
dial fluctuations is then dependent on angular resolution. 
Furthermore, we can develop rules of thumb for the in- 
strumental specifications required to unambiguously re- 
solve features such as the FoG effect and the bubble scale 
as we shall now see. 

4.3. Resolving redshift space distortions 

For FoG redshift space distortions, the critical scale 
where all fluctuations will be suppressed is associated 
with the pairwise velocity dispersion ap since that affects 
the small scale density contributions. For the ionization 
fluctuations, erf, typically places the suppression scale be- 
low that of the bubble scale and hence is less relevant. 
The pairwise velocity dispersion scale is independent of 
the bubble model and so leads to robust criteria for its 
resolution. 

Combining Figure [7| with Figure |H1 if we have both 
high frequency resolution Aiy < 0.03MHz and high an- 
gular resolution < 30", which means both suppress 
the power at smaller scales than the FoG effect does, we 
would be able to unambiguously detect the FoG effect 
in the flducial model with instrumental sensitivity that 
can extract mK level signals. More generally, we require 
that SD < ADo-p and DA^A9/y/8\n2 < AD^^ to ob- 
serve nonlinear redshift space distortions at small scales. 

Combining equation I^J, equation (|41|l . and equa- 
tion (|58|l . we obtain the rule of thumb on resolution re- 
quirements for observing FoG effect in 21cm experiments 



Af < Aiy„ = 0.034MHz 



l + z 
10 



30km/s 



(63) 



AO < A6I. = 8.9" 



l + z 
10 



1/2 



- 0.33 



V30km/s^ 
(64) 

neglecting the small cosmological dependence and assum- 
ing z 1. Note that redshift space distortions can only 
be resolved when both these criteria are satisfied. Like- 
wise for foreground removal, the FoG suppression be- 
comes the limiting factor only if both exceed this criteria. 

The characteristic scales defined above decrease with 
redshift. At 2; = 8 where we do all our calcu- 
lations, assuming ap — 30km/s, equation 163|l and 
equation (|64|l suggest the characteristic resolutions are 
Aiy^^ = 0.04MHz and AO^^ = 14". These resuhs are 
consistent with Figure [7| and Figure |H1 

4.4. Resolving HII regions 

We can also establish rulcs-of-thumb for the resolution 
of the HH regions or in our model the bubble size R. 



These supplement the considerations in Subsection 13.21 
for separability from the FoG scale. We will consequently 
consider the R = 10/i~^Mpc case where the two scales 
are themselves well separated. 
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Fig. 9. — Redshift space 21cm ID power spectra for different 
beam sizes and channel widths, assuming fiducial HII region size 
R = 10/i~^Mpc. For comparison, we also plot power spectra with 
perfect angular and frequency resolution (thick solid lines). 

Figure El shows how the redshift space power spectrum 
changes with different experimental resolutions. With a 
perfect instrument {A9 = and Av = 0) both the larger 
R scale and the FoG scale are resolved. However, as ei- 
ther increases, the small scale feature is replaced by the 
instrumental resolution. We can tell from Figure [HI that 
for Aiy < O.lMHz and A9 < 100", such large HH struc- 
tures can be distinguished. It is also interesting to see 
that when fixing Ad = 10", it does not help much to im- 
prove the frequency resolution beyond O.OlMHz, which 
implies the requirements for beam size and channel width 
are correlated in resolving structures, as we pointed out 
earlier in Subsection l4.2l 

More generally, we need both 6D <C R and 



D±A9/V8]ii2 < i? to resolve HII structures, 

10 ) \Mpc) 



Av < AvH = 0.14MHz 



(65) 



AO < AOr = 36.4" 



1 - 0.33 



10 

l + z 



1/2 



VMpcy 

(66) 

Either resolution not satisfying the above conditions will 
result in a smearing scale larger than HII scale so that 
the feature will be obscured. 

Furthermore, because of the existence of redshift space 
distortions, R and Da can lead to degenerate effects as 
we have already discussed in Subsection 13.21 Thus the 
required resolutions for observing any structures, either 
HII patches, or finger-of-God features, are the results of 
minimizing equation H63() and equation (|65|l . plus equa- 
tion 164|1 and equation H66|) . respectively. 
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5. CONSIDERATIONS FOR EXPERIMENTS 
5.1. An ultimate experiment 

In order to fully take advantage of all data, and go 
beyond the foregr ound d ominated region to remove fore- 
ground efficiently l|Wan g et al. 2005) , we would of course 
like to get as good angular and frequency resolution as 
possible. On the other hand, with a reasonable ioniza- 
tion model like that established here there are guides 
to when an instrument will recover essentially all of the 
useful cosmological signal. 

In the absence of redshift space distortions, we have 
seen that the predicted temperature fluctuations con- 
tinue to rise in rms as the frequency resolution increases 
due to density fluctuations outside of the ionization bub- 
bles. This would imply that to be most immune to 
smooth foregrounds in frequency space, one should con- 
centrate the instrumental sensitivity in as small a fre- 
quency band and field as possible. However, we have seen 
that FoG smoothing in the radial direction will generi- 
cally introduce a cutoff in the observed power and set 
a lower limit of Ai/ ~ lOkHz on the required frequency 
resolution. To match this frequency resolution one would 
desire an angular resolution of A9 = 10" although there 
are gains all the way out to A9 = 1". 

TABLE 1 

Resolution requirements for future "ideal" experiments 



Resolve 


R > ADap 


R < ADcTp 


R 


Ae^p <Ae< A0fl 

Au^p < Au < AuR 


Not applicable 


FoG 


Ad < Ae^p{< ABr) 
Au < AUap{< Aur) 


Ae < Ae^p 

Au < AUap 


Both 


Ad < Ae^p{< ABr) 
Au < Au^p{< Aur) 


Not applicable 



The requirements on resolving ionization bubbles are 
more model dependent. The best case scenario is that the 
HII patches are larger than the characteristic scale intro- 
duced by the velocity dispersion as one would expect near 
the end of reionization. Under this scenario, HII patches 
can be resolved without ambiguities from redshift space 
distortions and the experimental requirements can be re- 
laxed to satisfy A6' < A^/? and Aj/ < Avu. In the un- 
fortunate scenario when R < ADo-p, £-9-, if reionization 
ends with the creation of many small ionization regions 
or the observation redshift is too high, radial resolution 
of HII structures will be difficult even with perfect ex- 
perimental resolutions. 

Table ^ summarizes our angular and frequency resolu- 
tions considerations. Note that the signal to noise will 
also have to be high enough to detect the ^10 mK sig- 
nals in question. 

5.2. First generation experiments 

The first generation 21 cm experiments will be far from 
the ultimate angular resolution at the required signal to 
noise levels considered in the previous section. Wc dis- 
cuss here what the impact of redshift space modeling will 
be for these experiments. 



The upcoming MWA (Mileura Widefield Array) ex- 
periment^ has 8kHz frequency resolution which is in fact 
sufficient for even our ultimate experiment. At 200MHz 
(z — 6.1), its beam size is 3.4'. The characteristic reso- 
lutions at this redshift for resolving HII structures and 
observing FoG effects are given by equations 1)63(1 - 1166(1 as 
AiyR = 0.17MHz(i?/Mpc), A0r = 60"(i?/Mpc), Aiy^^ = 
0.05MHz(crp/30kms-i), and AO^^ = 17"(crp/30kms-i). 
These results suggest that MWA will be able to observe 
HII patches larger than ^ 3Mpc. It also implies that the 
FoG redshift space distortions will not be the limiting 
factor in foreground removal as long as the average ve- 
locity dispersion at that time is less than 350km/s, as 
expected in any reasonable model. 

At lOOMHz (z = 13.2), MWA can achieve an an- 
gular resolution ~ 6.8'. The characteristic resolu- 
tions at this redshift are Avb. = 0.12MHz(i?/Mpc), 
AdR = 50"(i?/Mpc), Aiy^^ = 0.02MHz(crp/30kms-i), 
and Ae„^ = 10"(crp/30kms"i). Similarly, it wiU not 
be able to resolve HII regions if they are smaller than 
^ 8Mpc and FoG distortions are even less of a limiting 
factor for foreground removal. 

Another upcoming experiment LOFAR"*, has frequency 
resolution 4kHz, angular resolution 3.1' for virtual core 
setup at 200MHz, and 5.2' at 120MHz. At 200 MHz (z = 
6.1), it will be able to observe HII regions larger than 
^ 3Mpc, yet will not be affected by FoG distortions as 
long as CTp < 300km/s. Similarly, at 120MHz (z = 10.8), 
it will resolve HII regions which are larger than ^ 6Mpc. 

Improving the angular resolution with longer base- 
lines can only improve these numbers before the char- 
acteristic scale for FoG distortions is reached ADu = 
0.3Mpc((T/30kms~^). If HII bubbles are smaller than 
this scale at the corresponding redshift, we are unable to 
unambiguously resolve them from their radial structure 
no matter how good the resolutions are. This complica- 
tion may impact the next generation experiment beyond 
MWA and LOFAR. 

In conclusion, for both LOFAR and MWA, the com- 
plication arising from nonlinear redshift space distortions 
in foreground cleaning should be minimal. Both experi- 
ments can detect large HII structures and their correla- 
tions at the end of reionization epoch. Yet it could be 
difficult for them to distinguish smaller patches expected 
at earlier times. 

One should also bear in mind that the actual abil- 
ity of individual experiments to distinguish HII regions 
also depends on many other factors such as quality of 
foreground cleaning, treatment of systematics, signal-to- 
noise ratio, etc. Experimental resolution is only one of 
many factors that determine how much information from 
the reionization epoch can be extracted from the obser- 
vational data. 

6. SUMMARY 

In this paper, we construct a self-consistent analytic 
model for 21cm brightness temperature power spectrum 
in the observed redshift space domain. In its simplest 
form, our model has two input parameters at each red- 
shift, the average ionization fraction and HII bubble size. 

3 littp://web. haystack. mit.edu/MWA/MWA. html 
* http://www.lofar.org 



11 



It can be readily generalized to consider a distribution of 
bubble sizes with various profiles while maintaining the 
proper scalings with ionization fraction. Though still a 
toy model of reionization, it is easy to apply to rapidly 
explore alternatives and has the flexibility to accommo- 
date various physical reionization mechanisms, including 
those based on atomic line cooling. 

Utilizing this model, we study the ID power spectrum 
along the line of sight. Given that foreground contami- 
nation is expected to be smooth in frequency space, the 
ID structure of the brightness temperature is the most 
robust observable and a likely stepping stone for con- 
structing the foreground cleaned 3D maps. 

We show the existence of redshift space distortions will 
eventually limit our ability to observe epoch of reioniza- 
tion. When the average size of HII regions is smaller 
than the characteristic scale of the average velocity dis- 
persion, as is likely toward the beginning of reionization, 
they can not be radially resolved. Furthermore at the low 
ionization levels typical of the beginning of reionization, 
the brightness fluctuations will be dominated by density 
fluctuations. Features due to the presence of HII regions 
may be masked by even relatively small uncertainties in 
the redshift space distortions. 

Combining the radial structure induced by the HII 
bubbles and nonlinear redshift space distortions, we out- 



line criteria for planning the angular and frequency res- 
olution of experiments. Even in the ID spectra, angu- 
lar resolution enters by eliminating contributions from 
modes that are not purely radial thus degrading the sig- 
nal. It is a serious limiting factor for the first generation 
experiments such as MWA and LOFAR. For such ex- 
periments, only HII bubbles that are larger than a few 
Mpc can be radially resolved even with ideal frequency 
resolution. On the other hand, nonlinear redshift distor- 
tions will not seriously affect such experiments. They do 
however suggest that for future experiments a frequency 
resolution of ~ lOkHz and an angular resolution of 10" 
will be sufficient to extract most of the information from 
the radial power spectrum near the end of reionization. 
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APPENDIX 

DISTRIBUTION OF BUBBLE SIZES 

A distribution of bub ble sizes can be straightforwardly accommodated in the framework of our model following 
IMortonsoTi fc"Hul (|2006^. Let us denote the probability of obtaining a bubble of radius within dR of R as P(R)dR. 
The ionization fraction associated with this bubble distribution is given by 

Xe = l-e-"''<^''^ (Al) 

where the average bubble volume is 

(Vb) - J dRP{R)Vb{R) . (A2) 

The window functions in one and two bubble term power spectra calculations are also replaced by average values. For 
two bubble terms, Wr in equation becomes 

(Wb) = 7^ f dRP{R)Vb{R)WR . (A3) 



{Vb 

For one bubble terms, in equation H3U|) is replaced by 

(Wr) = ^ / dRP{R)Vb^iRWl ■ (A4) 
Note that the approximation to the convolution term in equation (|34ll still holds with cr^ becoming 



(^Jj) = y dRP{R)Vb'iR)aj, . (A5) 

There are two qualitative effects of including a distribution. The averaging process removes the ringing of the 3D 
power spectra and broadens the feature associated with the bubble size. It also can change the relative strengths of 
the one and two bubble contributions. For example, since 

1™ {w^y ^ (u,)2 ' ^^^> 

the larger bubbles in the distribution enhance the power in the one bubble term at low k. 
To make these considerations concrete, let us consider a lognormal distribution 



(lni?-lni?o + 3<;j/2)2 



2<R 



(A7) 
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where the offset is defined so that (14) = Vq = 47r i?n/3 and similarly, (Vh) = Vnex-p(9a^„ p) . Numerical simulations 
suggest tha t a re asonable range for a\nR is [0.5, 1] ijFurlanetto. McQuinn. fc Hernauist II200?tD . 

In Figure 101 we show the ID and 3D real space power spectra for several values of din r and two choices of Rq. In 
the 3D case, all oscillations are averaged away and also the bubble transition takes place over a wide range of scales 
extending out to Ini? « lni?o + 3o'fj^^/2 due to the scaling of {Vff). Analogous effects occur for the ID power spectra 
and the shift in the transition scale toward larger scales can make the bubble transition slightly more prominent. 
Nevertheless the qualitative conclusions of the main paper remain unchanged once an allowance has been made for 
the effective characteristic bubble scale. 
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Fig. AlO. — Dependence of 3D (top) and ID (bottom) real space power spectra on the width of the bubble radii distribution <y\-^ji for 
two choices of i?o = 0.3Mpc//i (left) and 10Mpc//i (right). Note that oscillations in the 3D spectra are averaged out and the feature due 
to the bubbles is drawn out and extends to larger scales. Models are evaluated at 2; = 8 with = 0.8. 
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